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Abstract. The CRESST-II cryogenic Dark Matter search, aiming at detection of WIMPs via elastic scat- 
tering off nuclei in CaW04 crystals, completed 730 kg days of data taking in 2011. We present the data 
collected with eight detector modules, each with a two-channel readout; one for a phonon signal and the 
other for coincidently produced scintillation light. The former provides a precise measure of the energy 
deposited by an interaction, and the ratio of scintillation light to deposited energy can be used to dis- 
criminate different types of interacting particles and thus to distinguish possible signal events from the 
dominant backgrounds. 

Sixty-seven events are found in the acceptance region where a WIMP signal in the form of low energy 
nuclear recoils would be expected. We estimate background contributions to this observation from four 

sources: 1) "leakage" from the e/7-band 2) "leakage" from the a-particlc band 3) neutrons and 4) ^''^Pb 
recoils from ^^'^Po decay. Using a maximum likelihood analysis, we find, at a high statistical significance, 
that these sources alone are not sufficient to explain the data. The addition of a signal due to scattering 
of relatively light WIMPs could account for this discrepancy, and we determine the associated WIMP 
parameters. 

PACS. 95.35.-|-d Dark Matter, WIMP - 07.20.Mc Low-temperature detectors - 29.40.Mc Scintillation 
detectors, CaW04 - 29.40.Vj 



1 Introduction 

The nature of Dark Matter remains one of the outstand- 
ing questions of present-day physics. There is convincing 
evidence for its existence on all astrophysical scales and 
many theories predict particle candidates that may be able 
to explain its composition. However, in spite of numerous 
attempts, Dark Matter particles have not been unambigu- 
ously detected so far. 

Several experiments currently aim for direct detection 
of Dark Matter, mostly focusing on a particular class of 
particles, the so-called WIMPs (Weakly Interacting Mas- 
sive Particles). WIMPs nowadays are among the most in- 
vestigated and best motivated candidates to explain Dark 
Matter. If they exist, they could be present in our galaxy in 
the form of a halo, constituting the majority of the galac- 
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tic mass. Rare interactions with ordinary matter would 
then possibly be detectable in earthbound experiments. 

One such project is CRESST-II (Cryogenic Rare Event 
Search with Superconducting Thermometers), located at 
the Laboratori Nazionali del Gran Sasso in Italy. In this 
experiment we aim for detection of the WIMPs via their 
scattering off nuclei. The main challenges of this kind 
of measurement are to detect the tiny amounts of re- 
coil energy transferred to the nucleus (0(10 keV)), and to 
achieve sufficient background suppression to be sensitive 
to the extremely low rate of anticipated interactions (not 
more than a few tens of events per kilogram and year). 

To meet these requirements, CRESST uses cryogenic 
detectors in a low-background environment. Furthermore, 
the scintillating target material CaW04 allows for a dis- 
crimination of the type of interacting particle. In this way, 
potential rare WIMP interactions can be distinguished 
from events which were induced by the dominant radioac- 
tive backgrounds. 
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In this article, we report on the latest data obtained 
with the CRESST setup. They comprise a total net ex- 
posure of 730 kg d, collected between 2009 and 2011. The 
structure of the paper is as follows: Section [2] starts with 
a brief introduction to the experimental setup. Section [3] 
gives the new aspects of the present run and summarizes 
the observations. The choice of an appropriate acceptance 
region and the possible backgrounds relevant in this region 
are discussed. In Section [ij we give a qualitative descrip- 
tion and estimation of the backgrounds. Section [5] then 
describes the quantitative treatment in the framework of 
a maximum likelihood analysis. The results of this analy- 
sis are then discussed in Section [6j Finally, Section [7] gives 
the outlook for future runs of the experiment. 

2 Experimental Setup 

A detailed description of the CRESST-II setup was the 
subject of earlier publications pi, f?. Here, we restrict our- 
selves to a few key aspects relevant for the discussion of 
the new results. 

2.1 Target Material 

As a target for WIMP interactions, CRESST uses scin- 
tillating CaW04 crystals. They have a cylindrical shape 
(40 mm in diameter and height) and weigh about 300 g 
each. The current experimental setup can accommodate 
up to 33 of these crystals, constituting a maximum target 
mass of about 10 kg. 

Under the usual assumption of coherent WIMP scat- 
terings off nuclei, the scattering cross section contains an 

enhancement (with the mass number A of the target 
nucleus). In this case one expects that the total scatter- 
ing rate in CaW04 is dominated by interactions with the 
heavy tungsten nuclei. However, due to kinematics a heav- 
ier nucleus tends to receive a smaller recoil energy, and in 
a detector with a finite energy threshold the other con- 
stituents can also become relevant despite the coherence 
effect. To illustrate this point. Fig. [l] shows, as a function 
of the WIMP mass, the contributions of the three ele- 
ments in CaW04 to the total rate of WIMP interactions 
in the energy interval 12 to 40keV, the range typical of 
the CRESST detectors. 

For low WIMP masses, up to about 12GeV, the scat- 
terings off tungsten are completely below threshold, and 
oxygen and calcium recoils give the only possible Dark 
Matter signal. On the other hand, above WIMP masses 
of about 30 GeV, tungsten completely dominates. When 
looking for possible WIMP interactions we therefore con- 
sider recoils of all three types of nuclei in CaW04, in or- 
der to be sensitive to the largest possible range of WIMP 
masses. 

2.2 Phonon and Light Detectors 

In order to be able to detect the low energy nuclear recoils 
the target crystals are operated as cryogenic calorimeters 




WIMP mass [GeV] 



Fig. 1. Contributions of the three types of nuclei present in 
a CaW04 target to the total rate of WIMP interactions, as 
a function of the WIMP mass and for a cross-section of 1 pb, 
assuming coherent ~ interactions. The calculation assumes 
a detector with a sensitive energy range between 12 and 40 keV. 



at temperatures of about 10 mK. The energy deposited by 
an interacting particle is mainly converted into phonons, 
which are then detected with a transition edge sensor 
(TES) . We thus denote the target crystals with their TES 
as phonon detectors. 

The TES is a thin tungsten film evaporated onto the 
crystal, with the temperature stabilized in the transition 
from the normal to the superconducting state. The tiny 
change of the film temperature (O(iiK)) induced by the 
absorption of the phonons leads to a measurable change in 
resistance. This signal is read out by SQUID-based elec- 
tronics. The amplitude of the signal is a precise measure 
of the deposited energy. After the interaction, the crystal 
temperature relaxes back to the equilibrium state via a 
weak thermal coupling to the heat bath. 

In addition to the phonon signal, a small fraction of 
the energy deposited in the target crystal is converted into 
scintillation light. Each crystal is paired with a separate 
cryogenic light detector in order to detect this light sig- 
nal. Most of the light detectors are made from a silicon- 
on-sapphire wafer (a sapphire wafer of 40 mm diameter 
and 0.46 mm thickness, with a Ipm silicon layer on one 
side, which acts as photon absorber). As an alternative, 
some light detectors consist of pure silicon wafers of the 
same size. Similar to the target crystals, each light de- 
tector is equipped with a tungsten transition edge sensor 
to determine the energy deposited by the absorption of 
scintillation photons. 

A crystal and the corresponding light detector form 
a so-called detector module as shown in Fig. [2j Both de- 
tectors of a module are held by thin, silver-coated bronze 
clamps and are enclosed in a common, highly light reflec- 
tive housing in order to collect as much scintillation light 
as possible. The reflector is a polymeric foil which also 
scintillates. This will be discussed in detail in Section [2^41 
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Fig. 2. Schematic drawing of a CRESST detector module, con- 
sisting of the target crystal and an independent light detector. 
Both are read out by transition edge sensors (TES) and are 
enclosed in a common reflective and scintillating housing. 



2.3 Quenching Factors and Background Discrimination 

For each particle interaction, a detector module yields two 
coincident signals (one from the phonon and one from the 
light detector) . While the phonon channel provides a sen- 
sitive measurement of the total energy deposition in the 
target (approximately independent of the type of interact- 
ing particle), the light signal can be used to discriminate 
different types of interactions. To this end, we define the 
light yield of an event as the ratio of energy measured 
with the light detector divided by the energy measured 
with the phonon detector. We normalize the energy scale 
of the light channel such that 122 keV 7's from a ^^Co 
calibration source have a light yield of unity. With this 
normalization electron recoils induced either by /3 sources 
or by gamma interactions, generally have a light yield of 
about 1. On the other hand, a's and nuclear recoils are 
found to have a lower light yield. We quantify this reduc- 
tion by assigning a quenching factor (QF) to each type of 
interaction. The QF describes the light output expressed 
as a percentage of the light output for a 7 of the same 
deposited energy. 

Some quenching factors can be directly determined 
from CRESST data. For example, neutrons detectably 
scatter mainly off the oxygen nuclei in CaW04. The QF 
for oxygen can thus be determined from a neutron calibra- 
tion run which took place during the data taking discussed 
here. The result is 

QFo = (10.4 ±0.5)%. 

Moreover, the quenching factor for low energy a's can be 
found to be about 22 %, using a-events in the current data 
set. Similarly, the value for lead can be inferred to be 
around 1.4 %. Both measurements will be discussed below. 

Other types of interactions (in particular calcium and 
tungsten recoils) are not observed with sufficient statistics 
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Fig. 3. Illustration of background events due to surface con- 
taminations with ^^''Po. 



in CRESST, and their quenching factors must be deter- 
mined in dedicated experiments [3]: 
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Corresponding to these different values, there will be char- 
acteristic "bands" for the different particles or recoils in 
the light yield-energy plane. This allows for an excellent 
discrimination between potential signal events (expected 
to show up as nuclear recoils) and the dominant radioac- 
tive backgrounds (mainly e/7-events). 

Furthermore, it is even partially possible to determine 
which type of nucleus is recoiling. Such a discrimination 
is possible to the extent to which the different nuclear 
recoil bands in the light yield-energy plane can be sepa- 
rated within the resolution of the light channel. This then 
allows a study of potential WIMP interactions with dif- 
ferent target nuclei, in parallel in the same setup. Such a 
possibility can be particularly relevant for the verification 
of a positive WIMP signal, and is a distinctive feature of 
CRESST. 



2.4 Scintillating Housing 

As mentioned above, the housing of the detector mod- 
ules consists mainly of a reflecting and scintillating poly- 
meric foil. Making all surfaces in the vicinity of the de- 
tectors scintillating is important in discriminating back- 
ground events due to contamination of surfaces with a- 
emitters. The basic mechanism of this background is il- 
lustrated in Fig. |3] 

The most important isotope in this context is ^^°Po, 
a decay product of the gas ^^^Rn. It can be present on or 
implanted in the surfaces of the detectors and surround- 
ing material. The ^^°Po nuclei decay to ^"^Pb, giving a 
5.3 MeV a-particle and a 103 keV recoiling lead nucleus. 
It can happen that the lead nucleus hits the target crystal 
and deposits its energy there, while the a-particle escapes. 
Due to its low quenching factor, the lead nucleus can often 
not be distinguished from a tungsten recoil and thus can 
mimic a WIMP interaction. 

However, if the polonium mother nucleus was located 
on the surface of the target crystal or implanted in it (the 
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upper case in Fig. IS), the full 103 keV of the daughter 
nucleus plus a possible contribution from the escaping a- 
particle will be deposited in the target and the event will 
lie above the energy range relevant for the WIMP search. 

Another situation arises when the polonium atom was 
implanted in a surrounding surface. Then the daughter Pb 
nucleus can lose part of its energy on the way to the target 
crystal and appear in the energy range of interest (the 
lower case in Fig. [s]) . This possibly dangerous background 
can be rendered harmless if the surrounding surfaces are 
scintillating; in this case the escaping a-particle produces 
additional scintillation light when hitting those surfaces 
and the event will appear as high-light event in distinction 
to the low-light nuclear recoils. 

Hence the scintillation of the complete surroundings of 
the target crystals plays an important role. With the scin- 
tillating foil used as a module housing, currently the only 
non-scintillating surfaces inside the detector modules are 
the small clamps which hold the target crystals. In earlier 
runs, attempts were made to cover these clamps with scin- 
tillating layers as well, but these layers appeared to give 
rise to thermal relaxation events. The current module de- 
sign therefore avoids any scintillating (plastic) material in 
direct contact with the crystals. The price for this mea- 
sure, however, is the presence of several Pb recoil events 
with energies of 103 keV and below in the data set, as ex- 
pected from the above discussion. This background must 
therefore be taken into account in our analysis. 



2.5 Shielding 

Fig. |4] shows a schematic drawing of the whole CRESST 
setup, with the detector modules in the very center. The 
low temperatures are provided by a ^He-'*He dilution re- 
frigerator and transferred to the detectors via a 1.3 m long 
copper cold finger. The detector volume is surrounded by 
several layers of shielding against the main types of back- 
ground radiation: layers of highly pure copper and lead 
shield against 7-rays, while polyethylene serves as a mod- 
erator for neutrons. The inner layers of shielding are con- 
tained in a gas tight box to prevent radon from penetrat- 
ing them. In addition, an active muon veto using plastic 
scintillator panels is installed to tag muons. The veto sur- 
rounds the lead and copper shielding and covers 98.7% of 
the solid angle around the detectors, a small hole on top 
is necessary to leave space for the cryostat. 



2.6 Data Analysis 

We apply just a few basic quality cuts to the raw data 
in order to ensure that only valid events, with well- 
reconstructed energies in the phonon and light channel, 
are considered for further analysis. In particular, we re- 
quire that both the phonon and light detector of a given 
module were fully operational and running stably at their 
respective operating points at the time of an event. Data 
acquisition, readout, and the procedures for monitoring 
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Fig. 4. Schematic drawing of the CRESST setup. A cold fin- 
ger (CF) links the cryostat (CR) to the experimental volume, 
where the detectors are arranged in a common support struc- 
ture, the so-called carousel (CA). This volume is surrounded by 
layers of shielding from copper (CU), lead (PB), and polyethy- 
lene (PE). The copper and lead shieldings are additionally en- 
closed in a radon box (RB). An active muon veto (MV) tags 
events which are induced by cosmic radiation. 



detector stability, trigger efficiency as well as reconstruct- 
ing the deposited energy from the measured pulses are 
described elsewhere [T. 

For the final data set, we reject events coincident with 
a signal in the muon veto as well as those events with 
coincident signals in more than one detector module, since 
multiple scatterings are excluded for WIMPs in view of 
their rare interactions. 

An important aspect of the analysis concerns the 
bands in the light yield-energy plane. The e/7-band is 
highly populated due to the relatively high rate of com- 
mon backgrounds. This allows us to extract the position 
and, in particular, the energy-dependent width of this 
band directly from the measured data. The observed 
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Fig. 5. Data obtained with one detector module in a cali- 
bration measurement with an AmBe neutron source, with the 
source placed outside the lead shielding. The solid red lines 
mark the boundary of the calculated oxygen recoil band (10 % 
of events are expected above the upper and 10 % below the 
lower boundary). The vertical dashed lines indicate the lower 
and upper energy bounds of the WIMP acceptance region as 
will be introduced in Section [S] 

width is dominated by the light channel resolution com- 
pared to which the resolution of the phonon channel is 
much superior. This is understandable in view of the small 
fraction of the deposited energy appearing as light. 

We extract the resolution of the light channel as a func- 
tion of light energy by fitting the e/7-band with a Gaus- 
sian of energy dependent center and width. We note that, 
although the production of scintillation light is governed 
by Poisson statistics, the Gaussian model is a very good 
approximation in our regions of interest. This is because 
the e/7-events produce a sufficiently large number of pho- 
tons for the Poisson distribution to be well approximated 
by a Gaussian. On the other hand, for the quenched bands 
with low light yields, the Gaussian baseline noise of the 
light detector determines the resolution. 

The position and width of the bands other than the 
e/7-band can be calculated based on the known quenching 
factors discussed above and using the light channel reso- 
lutions obtained from the fit to the e/7-band. In order to 
get the width of a quenched band at a certain energy the 
light channel resolution for the actual light energy is used. 

To validate this calculation for quenched bands, we use 
the data from a calibration measurement with an AmBe 
neutron source placed outside the Pb/Cu shielding. We 
expect the neutrons to mainly induce oxygen nuclear recoil 
events. Fig. |5] shows the data obtained by one detector 
module in this measurement, together with the calculated 
central 80% band for oxygen recoils (10% of the events 
are expected above the upper and 10% below the lower 
boundary). 

Nuclear recoil events up to energies of about 300 keV 
are observed, with the spectrum falling off quickly towards 
high energies. In neutron-nucleus elastic scattering the re- 
coil energy of the nucleus is inversely proportional to its 
mass. Thus the highest energy recoils must be oxygen nu- 



clei. From the ratio of the mass numbers we then expect 
the highest energy of calcium recoils to be around 100 keV. 
Above 100 keV, we therefore have purely oxygen recoils, 
and the distribution fits well into the calculated oxygen 
band. Towards lower energies, the observed events are still 
in agreement with the prediction, although an increasing 
contribution from calcium recoils slightly shifts the center 
of the observed event distribution to lower light yields. 



3 The Latest Experimental Run 

3.1 Data Set 

The latest run of CRESST took place between June 2009 
and April 2011. It included a neutron test and 7-calibra- 
tions with ^''Co and ^"^^Th sources. In total, 18 detector 
modules were installed in the cryostat, out of which ten 
were fully operated. The remaining modules cannot be 
employed for a Dark Matter analysis, principally due to 
difficulties in cooling the light detectors. However, seven 
additional individual detectors (six phonon and one light 
detector) were still operated in order to tag coincident 
events (with signals in more than one module). 

One of the ten operational modules was equipped with 
a test ZnW04 crystal and we do not include it in this anal- 
ysis because of uncertainties in the quenching factors in 
this material. Another operational detector module had 
unusually poor energy resolution, with practically no sen- 
sitivity in the WIMP signal region, and was therefore ex- 
cluded from the analysis. The data discussed in this paper 
were thus collected by eight detector modules, between 
July 2009 and March 2011. They correspond to a total 
net exposure (after cuts) of 730 kg days. 

3.2 Observed Event Classes 

Fig. [6] shows an example of the data obtained by one de- 
tector module, presented in the light yield-energy plane. 

The e/7-events are observed around a light yield of 1. 
The calculated bands for a's, oxygen recoils, and tungsten 
recoils are shown The spread of a band at each energy is 
chosen so that it contains 80 % of the events, that is 10 % 
of the events are expected above the upper boundary and 
10 % of the events are expected below the lower boundary. 
This convention will be used throughout the following dis- 
cussion whenever we refer to events being inside or outside 
of a band. 

Beside the dominant e/7-background, we identify sev- 
eral other classes of events: 

Firstly, we observe low energy a's with energies of 
100 keV and less. They can be understood as a conse- 
quence of an a-contamination in the non-scintillating 
clamps holding the crystals. If the a-particle has lost 

^ The calcium band is not shown for clarity. It is located 
roughly in the middle between the oxygen and the tungsten 
bands. 
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Fig. 6. (Color online) The data of one detector module (Ch20), 
shown in the light yield vs. recoil energy plane. The large num- 
ber of events in the band around a light yield of 1 is due to 
electron and gamma background events. The shaded areas in- 
dicate the bands, where alpha (yellow), oxygen (violet), and 
tungsten (gray) recoil events are expected. Additionally high- 
lighted are the acceptance region used in this work (orange), 
the reference region in the a-band (blue), as well as the events 
observed in these two regions. See text for discussion. 



most of its energy in the clamp before reaching the target 
crystal, it can appear at low energy. The rate of such 
a-events differs by some factor of two among the detector 
modules (see Section |42]) . 

Secondly, Fig. [6] shows a characteristic event popula- 
tion in and below the tungsten band around 100 keV. This 
is present in all detector modules, albeit the number of 
events varies. This population can be attributed to the 
lead nuclei from ^^°Po a-decays on the holding clamps 



(see Section 2.4 1. The distribution of these events exhibits 
a low-energy tail, with decreasing density towards lower 
energies. In spite of this decrease, there are detector mod- 
ules (the ones with a high population of such lead events) 
in which the tail visibly reaches down to energies as low 
as a few tens of keV. 

Finally, low energy events are present in the oxygen, 
(calcium,) and tungsten bands at energies up to a few tens 
of keV, i.e. in the region of interest for the WIMP search. 
These events will be the main focus of our discussion in 
the following. We start by defining the acceptance region 
on which the discussion will be based. 



3.3 Acceptance Region 

Depending on the mass of a possible WIMP, any of the 
nuclei in CaW04 can be a relevant target for WIMP scat- 
tering as discussed above. We therefore choose our accep- 
tance region such that it includes all three kinds of nu- 
clear recoils: it is located between the upper boundary of 
the oxygen band and the lower boundary of the tungsten 
band. This selection automatically includes the calcium 
band. 



modulo 


Sr,'" [keV] 


acc. events 


Ch05 


12.3 


11 


Ch20 


12.9 


6 


Ch29 


12.1 


17 


Ch33 


15.0 


6 


Ch43 


15.5 


9 


Ch45 


16.2 


4 


Ch47 


19.0 


5 


Ch51 


10.2 


9 


total 




67 



Table 1. Lower energy limits E^^^ of the acceptance regions 
and the number of observed events in the acceptance region of 
each detector module. 



We restrict the accepted recoil energies to below 
40keV, since as a result of the incoming WIMP veloc- 
ities and nuclear form factors, no significant WIMP signal 
is expected at higher energies. On the other hand, to- 
wards low energies the finite detector resolution leads to 
an increasing leakage of e/7-events into the nuclear recoil 
bands. We limit this background by imposing a lower 
energy bound i?acc" each detector module, chosen such 
that the expected e/7- leakage into the acceptance region 
of this module is one event in the whole data set. Due to 
the different resolutions and levels of e/7-background in 
the crystals, each module is characterized by an individ- 
ual value of Table 111 lists the values of E"^!" for all 
modules. 

An example of the resulting acceptance region is shown 
(orange) in Fig. [6] and the events observed therein are 
highlighted. In the sum over all eight detector modules, 
we then find 67 accepted events, the origin of which we 
will discuss in the following. Table [T] shows the distribu- 
tion of these events among the different detector modules. 
Since fi'acc' "^^11 as the width of the bands are module- 
dependent, different modules have different sized accep- 
tance regions and thus different expectations with respect 
to background and signal contributions. 



3.4 Backgrounds in the Acceptance Region 

With the above choice of the acceptance region, four 
sources of background events can be identified: 

1. leakage of e/7-events at low energies, 

2. a-events due to overlap with the a-band, 

3. neutron scatterings which mainly induce oxygen recoils 
in the energy range of interest, and 

4. lead recoils from a-decays at the surface of the clamps, 
degraded to low energy. 

In the following, we estimate the contribution of each of 
these backgrounds and finally investigate a possible excess 
above this expectation. When present, such an excess may 
be the result of WIMP scatterings in our detectors, or of 
course an unsuspected background. 
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We estimate the backgrounds and the contribution of 
a possible WIMP signal in the framework of a maximum 
likelihood fit that allows a treatment with all relevant pa- 
rameters and their uncertainties simultaneously. However, 
before introducing this rather abstract formalism in Sec- 
tion[5] we give a qualitative discussion of the backgrounds, 
with the aim of clarifying the basic arguments and as- 
sumptions. 



4 Qualitative Background Discussion 

4.1 e/7-Background 

The lower energy bound of the acceptance region is chosen 
such that we expect a leakage of one background e/7-event 
per detector module. These events are expected to appear 
mainly at the low energy boundary of the acceptance re- 
gion where the overlap between the bands is largest and 
towards the upper boundary, closer to the e/7-band. The 
quantitative modeling is discussed in Section [5. 2[ 

4.2 a-Background 

Since the a-band has some overlap with the acceptance 
region, some low energy a-events may be misidentified as 
oxygen or even calcium or tungsten recoils. This will lead 
to a certain expectation n"^^. of background events in the 
acceptance region of each module. 

The energy spectrum dNa/dE of the low energy events 
in the a-band appears to be flat within the available statis- 
tics. In order to estimate dNa/dE, we first select a refer- 
ence region of the a-band which is free of overlap to any 
other band. 

An example of such a reference region is highlighted 
(blue) in Fig. [6] To obtain reasonable statistics, we have 
chosen the relatively large energy range of 100 keV for 
the reference region. The low energy limit Ef^^ of the 
reference region is chosen as low as possible, while tak- 
ing into account the increasing e/7-leakage into the a- 
band towards low energies. Low i?™ are desirable, since 
the reference region should naturally be close to the ac- 
ceptance region to minimize extrapolation errors. On the 
other hand, influences of the e/7-background on the a- 
estimate should be minimized and the value of i?™™ was 
thus chosen such that only 0.1 counts of e/7- leakage are 
expected in the whole reference region for each detector 
module. Since the width of the bands varies from module 
to module, each module has its own value of These 
values are listed in Table [2] 

For our first rough estimate of the alpha background 
in the acceptance region, we simply count events in the 
reference region. We then assume constant dNa/dE and 
calculate the ratio of a-events expected in the acceptance 
region to those expected in the reference region. Scaling 
the observed number of events in the reference region by 
this ratio, we arrive at an estimate of the a-background 
in the acceptance region. 



module 


E^j" [keV] 


"-"of 


"-acc 


Ch05 


21.7 


17 


1.6 


Ch20 


21.3 


14 


1.5 


Ch29 


21.7 


14 


1.2 


Ch33 


28.3 


8 


0.9 


Ch43 


29.7 


8 


0.6 


Ch45 


24.7 


5 


0.8 


Ch47 


32.2 


9 


1.2 


Ch51 


18.3 


18 


1.4 


total 




93 


9.2 



Table 2. Lower energy limits -E™f" of the a-reference regions, 
observed alpha counts n'^^i in the reference regions, and the 
resulting (rough) estimates of the alpha background n^cc in 
the acceptance regions of all detector modules. 



Table [2] summarizes the observed alpha counts nf^f in 
the reference region and the resulting estimates of the 
alpha background n^^.^. in the acceptance region of each 
module. This results in a total expected a-background of 
about 9.2 events. 

In the likelihood analysis of Section [5j we will relax 
the assumption of constant dNa/dE and also allow for a 
linear term in the a-energy spectrum. We will, however, 
see that the simple estimate given here is quite close to 
the one obtained with the more sophisticated analysis. 



4.3 Neutron Background 

4.3.1 Introduction 

Throughout our discussion, we distinguish two different 
classes of neutron production mechanisms: 

Firstly, free neutrons can be emitted by radioactive 
processes, in particular spontaneous fission (s.f.) of heavy 
elements or (a, n) reactions on light nuclei. Such neutrons 
typically have energies up to a few MeV, for which the 
polyethylene shielding provides a very efflcient modera- 
tor. Moirte Carlo simulations suggest that neutrons from 
s.f. and (a, n) reactions in the rock outside the experiment 
thus only constitute a negligible background at the level 
of 10~^ events per kgd Therefore, such neutrons are a 
possibly relevant background only if they are emitted in- 
side the neutron shielding, e.g. by s.f. of ^ssy in the lead 
shielding, or by (a, n) reactions or s.f. in the copper shield- 
ing. We consider such neutrons here, even though Monte 
Carlo simulations predict only a negligible background in 
the acceptance region due to these processes, at a level of 
10~'^ events per kgd [5]. 

Secondly, neutrons can also be produced by muons, ei- 
ther in the lead or copper shielding or in the surrounding 
rock. In the former case, the muon will mostly be tagged 
by the muon veto enclosing the Pb/Cu shielding. How- 
ever, there is a small probability that the muon is missed 
by the veto because of the hole on top of the setup, as 
mentioned above (cf . Fig. W| . Such a muon may create a 
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shower of neutrons inside the PE shielding which then 
reach the detectors. On the other hand, muon-induced 
neutrons created outside the neutron shielding may pene- 
trate the polyethylene layers if they are energetic enough. 
Such high energy neutrons then have a high probability to 
scatter inelastically in the Pb/Cu shielding and to create 
secondary neutrons and gammas. Ultimately, this leads to 
events with similar characteristics as when the shower is 
directly induced by a muon inside the Pb/Cu shielding. 



4.3.2 Experimental Information 

It is a characteristic feature of neutrons that they can lead 
to coincident events in more than one detector module 
at the same time, with at least one module registering a 
nuclear recoil. A given neutron source will thereby lead 
to events with a characteristic ratio between single and 
coincident scatterings. If this ratio is known, the observed 
coincidences can be exploited to estimate the expected 
number of single neutron (background) events. This is the 
basic concept pursued in the following. 

We shall base our discussion on two sources of informa- 
tion. One is the results of the neutron test with an AmBe 
neutron source already mentioned above. The second is 
the examination of events in coincidence with incoming 
muons, i.e. muon veto triggers accompanied with a signal 
in the acceptance region of a module. These two "cali- 
bration measurements" are used to infer the properties of 
neutron backgrounds originating from ambient neutrons 
and from muons escaping the veto and interacting in the 
apparatus. We emphasize that Monte Carlo or other ex- 
ternal information is helpful, but does not enter into our 
quantitative estimates. 

For each type of neutron mechanism, we use the cor- 
responding calibration to infer the typical structure of the 
events with respect to their multiplicity, defined as the 
number of detector modules triggering at the same time 
(within a time window of 5 ms) . For muon-induced neu- 
trons, we observe 40 events coincident with a muon veto 
trigger and with at least one detector module having a 
signal in its acceptance region. The multiplicities of these 
events are shown in the top histogram of Fig. [7] On the 
other hand, the multiplicities of the events induced by the 
AmBe neutron source (placed at various positions outside 
the Pb/Cu shielding but inside the neutron shielding) are 
given in the bottom histogram of Fig. [7j 

Events due to muon-induced neutrons show an obvi- 
ously higher average multiplicity than events from the 
neutron test source. This evidently results from muon- 
induced cascades, leading to neutrons and gammas in dif- 
ferent detector modules at the same time. 

With this information at hand, we turn back to the 
Dark Matter data set. In addition to the 67 accepted 
events, these data contain three events in which several 
detector modules triggered in coincidence, with at least 
one module registering an event in its acceptance region. 
Two of these events have a multiplicity of three (i.e. three 
modules triggered), while in the third event five modules 
triggered in coincidence. Such coincidences which include 
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Fig. 7. Multiplicity (i.e. the number of modules which trig- 
gered in coincidence) of events which include a nuclear recoil 
in the acceptance region of at least one detector module. Top: 
muon-induced events, bottom: events from an AmBe neutron 
source. 



at least one nuclear recoil must involve a neutron, and 
perhaps 7's. Accidental coincidences may be neglected in 
view of the low overall counting rate. 

Based on these three coincident events we can use the 
calibration information given above to scale up to the 
number of expected single scatters for each type of source. 

For a first rough estimate of the level of neutron back- 
ground, we start with the three observed coincidences and 
neglect their precise multiplicities. From the histograms in 
Fig.[7| we find that muon-induced neutrons are character- 
ized by a ratio of single to coincident scatterings of about 
0.5, while the neutron source test gives a ratio of about 
3.8. Consequently, if we assume that only muon-induced 
neutrons are present in the experiment, this would lead to 
an estimated number of single scatterings of 3 • 0.5 = 1.5. 
On the other hand, in case the neutron background purely 
comes from a radioactive source, the same estimate gives 
a background expectation of about 3 • 3.8 = 11.4 single 
events. 

In reality, both types of neutrons may be present, and 
the above limiting cases show the extremes between which 
the expected neutron background can lie according to this 
estimate. Our likelihood analysis of Section [5] takes into 
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account the more detailed information of the individual 
event multiplicities in order to clarify the contributions of 
the two types of neutron sources to the total background. 
We will, however, see that the result is compatible with 
the simple estimates of the limiting cases given here. 

An independent aspect of the neutron background con- 
cerns the corresponding recoil energy spectrum. Within 
our narrow accepted energy range, the energy spectra 
induced by the two types of neutron events are found 
to be very similar, according to the calibration data 
discussed above. The spectrum can be parametrized by 
a simple exponential dNn/dE oc exp{—E/Edec)- We 
determine the parameter i^dec from a fit to the spec- 
trum obtained in the AmBe neutron calibration run. In 
the energy range between 12keV to 40keV we obtain 
£;dec = (23.54 ± 0.92) keV. 

This similarity in the spectra induced by neutrons from 
the two quite different sources (in agreement with Monte 
Carlo results [5]) indicates how the Pb/Cu shielding sur- 
rounding the detectors will moderate an incoming neu- 
tron flux regardless of its origin. The primary spectrum of 
the neutrons is washed out by inelastic scatterings in the 
shielding. This finding supports our use of the results of 
the neutron calibration to estimate the effects of a gen- 
eral neutron background. The only exception to this ar- 
gument might be a neutron-producing contamination in 
close vicinity of the detectors. In this case, we would ex- 
pect a recoil spectrum reaching to much higher energies 
and fewer singles for a given number of coincidences. In 
this case, the application of our above calibration results 
would lead to a conservative neutron background estimate. 



4.4 Lead Recoil Background 

To illustrate the lead recoil background from ^^'^Po decay, 
Fig. [8] displays the data set of a different detector mod- 
ule as in Fig. [6l Compared to Fig. [6j a more prominent 
population of Pb recoils below the tungsten band is 
visible, with a rather long tail extending down to the ac- 
ceptance region. Since the lead band and the acceptance 
region overlap considerably, a leakage of some ^"^Pb events 
into the acceptance region cannot be excluded. 

For an estimate of this background, we follow a sim- 
ilar strategy as for the a-background. We define a refer- 
ence region for each detector module which contains pre- 
dominantly ^°^Pb recoils, and model the spectral energy 
density dNp^/dE in this region. This model is then ex- 
trapolated into the energy range of the acceptance region. 

As a reference region, we choose the lead recoil band 
at energies above the acceptance region, where a possible 
WIMP signal cannot contribute. In some detector modules 
with wider bands, the lead band still overlaps with the 
oxygen band around the lower edge of this energy range. 
In this case, we additionally restrict the reference region 
to the lower part of the lead band without overlap with 
the oxygen band in order to be independent of possible 
neutron-induced events on oxygen. The event distribution 
of the Pb recoils peaks at the full lead recoil energy of 
103 keV and the upper boundary of the reference region 




60 80 
Energy [keV] 

Fig. 8. (Color online) The data of detector module Ch51, 
shown in the light yield vs. recoil energy plane. Again, the 
shaded areas indicate the bands, where alpha (yellow), oxygen 
(violet), and tungsten (gray) recoil events are expected. Ad- 
ditionally highlighted are the acceptance region (orange), the 
region where lead recoils with energies between 40 and 90 keV 
are expected (green), and the events observed in these regions. 
The highlighted lead recoil region (green) serves as a reference 
region for estimating the ^"^Pb recoil background. 



module 




Ch05 


17 


Ch20 


6 


Ch29 


14 


Ch33 


6 


Ch43 


12 


Ch45 


15 


Ch47 


7 


Ch51 


12 


total 


89 



Pb 
rcf 



Table 3. Observed counts n^^ in the lead reference regions of 
the detector modules. 



is set at 90keV so that it covers the low energy tail. An 
example of the resulting reference region is highlighted 
in green in Fig. [sl Table [3] summarizes the counts n. 
observed in the reference region of each detector module 
Fig.lOlpresents the energy spectrum of the events found 
in the "^^Ph reference regions of all detector modules, but 
includes also lead recoils with higher energies to illustrate 
the peak at the full nominal recoil energy of 103 keV. In 
the energy range of the reference region (below 90keV), 
the tail of the distribution can be modeled by an expo- 
nential decay on top of a constant contribution: 



dM 



Pb 



dE 



(E) = Apb 



Cpb + exp 



£;-90keV 

decay 



(1) 



For a first rough estimate of the recoil background, 
we simply fit such a function to the spectrum of Fig. [9] 
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70 80 
Energy [keV] 

Fig. 9. Energy spectrum of events in the ^°^Pb reference region 
of all detector modules. Events occurring in the lead band at 
energies above the upper energy limit of the reference region 
are included to show the peak at the nominal recoil energy 
of 103 keV. The curve (red) is a fit of the histogram with an 
exponential plus a constant term, in the energy range of the 
reference region. 



The red line shows the result of this fit with the param- 
eters Apb = 4.53counts/keV, Cpb = 0.13 and E^^^^^ = 
13.72 keV. This model then needs to be extrapolated into 
the energy range of the acceptance region. 

To check the validity of this extrapolation we used the 
SRIM package [6] to simulate the energy spectra expected 
for three different depth distributions of the ^^°Po mother 
nucleus. These distributions were: an exponential profile 
with 3 nm decay length peaking at the surface, a uniform 
distribution in the volume, and finally the depth distribu- 
tion resulting from implantation due to preceding alpha 
decays. For the latter case the implantation profile was 
also calculated with SRIM, assuming that ^^^Rn is first 
adsorbed on the surface of the clamps holding the crys- 
tals, followed by two subsequent alpha decays to ^^"Po. 
The results of the three simulations are shown in Fig. [TO] 

An important result of the simulation is that none of 
the calculated spectra of the Pb recoils rises significantly 
towards low energies within the range of our acceptance 
regions. The simulated spectra for the uniform as well as 
the implantation profile are rather flat between 40 and 
90 keV compared to the data in Fig. [9j However, the en- 
ergy spectrum from the distribution peaking at the surface 
has a tail similar to that observed in the data. The curve 
in Fig. 10 is the result of a fit of Eq. [l] to this spectrum 
in the energy range of the reference region between 40 
and 90keV. Its exponential decay of E^^^ — 13.6 keV 
agrees within errors with the value obtained from the fit 
of the data in Fig. [9] Below 40 keV the curve in Fig. [T0| 
shows the extrapolation of this function into the accepted 
energy range. The very good agreement of extrapolation 
and simulated data in the energy range of the acceptance 
region justifies the use of this type of extrapolation for our 
estimate of the Pb background. 
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Fig. 10. Energy spectrum of '^''''Fb recoils calculated with the 
SRIM package for three different depth profiles of the alpha 
emitting ^^"Po parent nuclei: (top) distribution peaking at the 
surface with a 3 nm exponential decay length, (middle) uniform 
distribution throughout the volume, and (bottom) a depth pro- 
file resulting from adsorption of ^^'^Rn at the surface of the 
clamps, with a subsequent implantation by the two alpha de- 
cays which follow in the -^^^U decay chain. The full (red) curve 
is the result of a fit of Eq.[T]to the top spectrum, in the energy 
range of the reference region from 40 to 90 keV. The extrap- 
olation of this fit into the acceptance region below 40keV is 
shown. 



For a rough first result, we take a typical energy range 
for the acceptance region of 12 to 40keV and estimate the 
number of Pb recoils in it. From the extrapolation of the 
fit function in Fig. [9] we calculate about 17 events of ^''^Pb 
background in the acceptance region. Of course, this sim- 
ple estimate neglects small differences in the overlap of 
the lead band with the acceptance region in different de- 
tector modules, as well as the acceptance reduction in the 
reference region due to partial overlap of lead and oxygen 
bands. Nevertheless, the result is very close to the final 
value that we will obtain from the full likelihood analysis, 
which performs the background estimate module- wise and 
hence takes such differences into account. 



5 Maximum Likelihood Analysis 



In the previous section, the qualitative principles of our 
background estimates were discussed. This section ex- 
plains how we formulate these concepts quantitatively. 
We choose the framework of a likelihood analysis to esti- 
mate the unknown parameters of our backgrounds and as 
well as the ones of a possible signal. This formalism also 
allows us to take into account and propagate the corre- 
sponding uncertainties. We first give a general overview of 
the formalism before focusing in detail on the treatment 
of the different backgrounds in the current measurement. 
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5.1 General Concepts 

The maximum likelihood fit is based on a parameterized 
model of the backgrounds and a possible signal. The pa- 
rameters are then varied to find the values for which the 
model most likely reproduces the data. 

As we will see, this procedure allows one to take into 
account the actual position of the events in the light yield- 
energy plane, and not merely their presence in a broad 
acceptance region. Also, measurements from reference re- 
gions, outside the acceptance regions, may be introduced 
to help determine some of the parameters. Furthermore 
one may, in the framework of the method, give a quan- 
titative estimate of how well a best set of parameters is 
determined. 

The maximum likelihood method itself however gives 
no direct indication of the quality of the resulting fit. To 
deal with this question, we shall attempt to give a charac- 
terization of the quality-of-fit using the so-called p-value. 



5.1.1 Likelihood Function 

The basic input to the analysis is a model for each source 
X that may contribute to our accepted events, be it back- 
ground or a possible signal. Such a model describes the 
expected distribution of events in the (E, y) plane (where 
E denotes the recoil energy and y the light yield). 

We formulate such a model in terms of a density 
px{E^y\'Px) in the {E,y) plane. The number of events 
expected from source X in any subspace of the (E, y) 
plane is then simply given by the integral of px over 
this region, px depends on a set of unknown parameters 
summarized in the vector Tpx- These are the parameters 
to be varied. 

In our case, each such density function is the product 
of two components: the expected recoil energy spectrum 
dNx/dE and a function that describes the expected event 
distribution in the light yield coordinate at each recoil 
energy E. 

As discussed above, we use a Gaussian distribution for 
the light yield, where the center of the distribution is given 
by the respective quenching factor and the width is essen- 
tially determined by the resolution of the light channel. 

Since different detector modules have different resolu- 
tions, the densities px need to be defined for each module 
individually. We therefore add an additional index d for 
"detector module" to each density. Assuming that we have 
considered all significant sources of events, the total den- 
sity for module d is 



/(i?,y|p) = p^ 



pi 



Pncutron 



PPb 



(2) 



where we have taken into account the four backgrounds 
discussed above, plus a possible WIMP signal "x"- The 
vector p summarizes all the unknown parameters. 

We now have observed a set of events in the acceptance 
region of each module, located at the positions {Ei,yi). 
The functions p'' can be fitted to the event distribution 
in the respective module and, from this the most likely 



values of the unknown parameters can, in principle, be 
determined. The normalizations of the p depend on the 
variable parameters and give the total number of expected 
events M via 



p\E,y\p) dEdy 



acc. region 



(3) 



where the integral runs over the acceptance region of mod- 
ule d and thus yields the total expected number of ac- 
cepted events in this module. 

The p'^ are thus no probabilty densities, but they can 
be used in the so-called extended maximum likelihood for- 
malism U to formulate the likelihood function which then 
takes the form 



'Cacc(p) = 



L d 



exphA/-(p)], (4) 



with the outer product running over the detector modules 
and the inner one over all accepted events in the respec- 
tive module. The exponential factor takes into account the 
normalization of the p^. 

We emphasize that, in Q, we directly evaluate p'^ at 
the position of each observed event and no binning is in- 
volved. By fitting the density functions to the observed 
event coordinates, one makes use of the full available in- 
formation from the two-channel measurement. Using the 
experimental data, the fit then finds the most likely val- 
ues of the parameters, including those for the backgrounds 
and a possible signal. 



5.1.2 Reference Regions 

In practice, some of the unknown parameters in p are not 
sufficiently constrained by the observation in the accep- 
tance region alone. In such cases, additional observations 
("reference measurements") can be exploited which are 
made either in other regions of the measurement param- 
eter space (for example in the a-reference region intro- 
duced in the previous section), or in completely different 
experiments or measurement channels. Each such refer- 
ence measurement yields a separate likelihood function 
which constrains a subset p' of p, and, assuming all mea- 
surements are independent, the total likelihood is simply 
given by the product 



Aot(p) = >Cacc(p) • 'Ci.cfl(pi) • 'Crcf2(P2) 



(5) 



reference measurements 



This concept can also be exploited in order to take 
into account the effects of uncertainties on otherwise fixed 
quantities, like e.g. the quenching factors of O, Ca and W, 
which enter the analysis and have been measured previ- 
ously. To this end, these quantities are added to p as free 
parameters and, at the same time, an additional likelihood 
term is included for each of them which models the uncer- 
tainty of the quantity in question (e.g. a Gaussian around 
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the best estimate, with the width given by the error on 
this estimate). 

Finally, maximizing the total likelihood function £tot(p) 
leads to estimates for all unknown parameters p. Most 
of the time, however, one is only interested in a subset 
of these values, the others just being nuisance parame- 
ters required to construct the likelihood. In our case, the 
only parameters of interest are the WIMP-nucleon cross 
section awN and possibly (if we find a clear signal) the 
WIMP mass m^. Our aim is to derive confidence inter- 
vals for these relevant quantities, taking into account the 
nuisance parameters and their uncertainties. 

5.1.3 WIMP Parameters 

In a first step, we will thereby only concentrate on ctwn 
and ask whether our measurement indicates awN > sig- 
nificantly. All other parameters shall be summarized in 
the vector p'. This question corresponds to a test of the 
null-hypothesis ctwn — with the model. 

As a convenient test statistic, we employ the likelihood 
ratio ^ 

^(o-wN = 0) := — — ^— . (6) 

Here, (Twn and p' in the denominator are the maximum 
likelihood estimators of ctwn and p', respectively, and p' is 
the conditional maximum likelihood estimator under the 
condition (Twn = 0. yl is hence a measure of how "signal- 
like" our observation is, in the sense that small values 
indicate a considerably better description of the data with 
an allowed WIMP contribution than without. 

According to Wilk's theorem [8j, if (Twn — were the 
true assumption, then the quantity 

g := -21nyl(crwN = 0) (7) 

would have a ^^-probability distribution for one degree of 
freedom in the limit of large statistics. We have verified 
with a Monte Carlo simulation that this approximation 
holds well for the case of our analysis. This allows for a 
simple calculation of the statistical significance S", with 
which we can reject the null hypothesis ct-wn — when 
having observed a certain value ^obs: 

S = V?obs- (8) 

Higher values of S obviously imply smaller probabilities 
to get data as "signal- like" as the observed ones, although 
there is no true signal present. 

In a second step, if ctwn > can be established with 
sufficient significance, the above approach can be gener- 
alized and can be treated as a second parameter of 
interest. The aim is then to calculate a confidence region 
in the (m^, ctwn) plane. Such confidence regions are given 
by the contours on which the likelihood has decreased by a 
certain factor 5C with respect to the maximum, provided 
that all parameters other than and ctwn are refitted in 
each point. The value of 6C that yields the desired confi- 
dence level can thereby be obtained from a x^-distribution 



for two degrees of freedom. In fact, the same procedure can 
be applied to calculate confidence regions for any combi- 
nation of rt > 1 parameters, using the ^^-distribution for n 
degrees of freedom. The one-dimensional case is commonly 
known from its implementation in the MINOS program of 
the MINUIT software package [S]. 

5.1.4 p-Value 

As mentioned above, the significance resulting from the 
above likelihood ratio test is, in itself, not a good measure 
of the quality of the fit to the data. High significances do 
not automatically imply a good agreement between the 
fitted model and the data, so that the quality-of-fit has to 
be determined in an independent step. One possibility of 
doing so is described in [TD] and we adopt an analogous 
procedure here: 

The {E, ?/)-plane (or a subspace of it) is divided into 
two-dimensional bins, separately for each detector module. 
In each bin i, the expected number of events can be 
calculated by integrating the density p'^. To go with this 
expectation, we have certain numbers of observed events 
in the bins. We can then determine the total probability 
of this observation by assuming an independent Poisson 
process with expectation in each bin i and multiplying 
all corresponding Poisson probabilities. 

Let our real observation have such a probability Pobs- 
Independently, we can then generate artificial data sets 
where, in each bin i, the generated number of events is 
Poisson distributed with the expectation e^. We employ a 
Markov Chain Monte Carlo for a fast generation of these 
data sets. Each generated set s can then again be charac- 
terized by a certain probability Ps. From the distribution 
of the Ps , we finally determine the probability to get a data 
set with probability less or equal to Pobs- This is called the 
p- value and can be used to describe the quality-of-fit. 

Often, as in our case, the model used to generate the 
artificial data sets is the result of a previous fit to the real 
observed data, where parameter values have been tuned 
for an optimal description of this observation. The p- value 
as obtained above needs to be corrected for this bias and a 
procedure of doing so is suggested in [lOj . We perform this 
correction, but we note that its influence is naturally small 
when the number of free fit parameters is small compared 
to the number of bins considered. 



Having now established the formal framework of our 
analysis, we will focus in the following on the concrete con- 
struction of the densities as well as the required reference 
measurements for each considered source of events. 



5.2 e/7- Background 

The distribution of the e/7-events in the {E,y) plane is 
observable directly in the Dark Matter data set. Due to 
the large number of such events, statistical fluctuations 
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are negligible and we use the observed distribution di- 
rectly to fix the densities p'^{E,y). The fit of the energy- 
dependent Gaussian width of the e/7-band as outlined 
above thereby gives the parametrization of the light yield 
distribution of the events. On the other hand, we use a his- 
togram of the observed recoil energies (bin width 0.1 keV) 
as the expected energy spectrum and interpolate it to ar- 
rive at a continuous function. In the energy range of in- 
terest, the statistical uncertainties of the bin contents of 
this histograms are below 10~^. The resulting densities pij 
can therefore be treated as fixed contributions to the total 
density in each module, with no free fit parameters. 



5.3 a-Background 

As discussed above, the observation of events in overlap- 
free regions of the a-band indicates an approximately flat 
energy spectrum of the a-background in our data. Never- 
theless, to account for a possible small systematic energy 
dependence, we model the energy spectrum in each detec- 
tor module d as 



dN: 



dE 



HE) 



(9) 



introducing the fit parameters c„ and TOq. c„ is module- 
specific and describes the level of a-activity observed in 
the respective module, while we assume the possible slope 
ma to be common to all modules. It reflects the properties 
of the alpha energy loss before interaction in the target. 

Obviously, the above fit parameters are only weakly 
constrained by the events observed in the acceptance re- 
gion, where other sources contribute as well. Therefore, we 
additionally consider the a-reference region for each de- 
tector module as introduced in Section |4] It is this region 
that mainly influences the estimates of the a-parameters. 
Technically, we simply extend the inner product in Eq. Q 
to also run over the events observed in the reference re- 
gion and modify the integral in Ec^. ([s]) accordingly. Since 
the contributions from other sources are small in the refer- 
ence regions, this modification mainly leads to additional 
constraints on the a-parameters. 

A second benefit of using these reference regions is that 
the low-energy quenching factor for a-particles can be left 
free in the fit. This is particularly important, since our 
observation of alphas in this data set is the best measure- 
ment available for this quantity. The quenching factor is 
automatically constrained by the events observed in the 
reference region and the uncertainties of this measurement 
are directly built into the likelihood. Of course, in order 
to define the reference region, a reasonable starting as- 
sumption for the quenching factor is required. We obtain 
it simply from the one-dimensional light yield distribution 
of the low-energy a-events. If this starting assumption is 
not too far from the true value, the likelihood fit can give 
a good estimate for the quenching factor, even though it is 
based on the events previously selected under the assump- 
tion of a slightly different starting value. We find that the 
result is indeed very robust against reasonable changes in 
the starting value. 



We also confirmed that the result is very robust against 
some reasonable changes of the energy range of the refer- 
ence regions. A selection of the same energy range for all 
detector modules of 35 to 135 keV, or selecting a 150 keV 
wide reference region with the lower energy limits of Ta- 
ble [2] did not change results in any relevant way. 



5.4 Neutron Background 

For the neutron background, we consider the two types 
of neutron creation mechanisms discussed in the previous 
section. We have seen that both, a radioactive neutron 
source and muon-induced neutrons can be described by 
the same exponential energy spectrum of the correspond- 
ing recoil events. Our model is hence 



d^n,^, , , , { E 

——{E) = An(nrad + nmuon ) ' exp - — — 



(10) 



with E'dcc = 23.54 keV as given in the previous section. 
The amplitude A„ is a function of the total expected num- 
ber of neutrons, i.e. of the sum Tiiad + "-muon of accepted 
single neutron scatterings due to radioactive sources and 
muon-induced neutrons, respectively. Both these numbers 
are unknown parameters. These parameters shall here be 
defined as the sum of expected events in all detector 
modules. Nevertheless, ( jTo] ) describes the expected energy 
spectrum of neutron-induced events in each module. 

Although neutrons will scatter off all three types of tar- 
get nuclei in CaW04, Monte Carlo simulations show that 
more than 90% of all scatterings detected in the energy 
range of interest happen off oxygen. In our background 
model, we hence treat neutron events as oxygen recoils 
and use the corresponding quenching factor to describe 
their light yield distribution. The uncertainty of the oxy- 
gen quenching factor as given in Section |2.3| is included in 
the likelihood, assuming a Gaussian error. 

With also other sources of events contributing to our 
acceptance regions in the oxygen band, the unknown pa- 
rameters rirad and nmuon of the neutron background need 
to be constrained by an additional reference measurement. 
As outlined above, we exploit the ability to cause coin- 
cident events in more than one module to estimate the 
contribution from neutrons. 

Each of the two neutron creation mechanisms consid- 
ered here has a certain expected ratio of coincident events 
with multiplicity m > 1 to single scatterings: 



rad / muon 



exp. events with mult, m in acc. region 
exp. single events in acc. region 



The multiplicity histograms of Fig. [7] represent our mea- 
surement of these ratios. 

On the other hand, we have observed coincident events 
in our background data as outlined above and we will use 
them to derive constraints on the expected number of sin- 
gle neutron scatterings. To this end, we treat the occur- 
rence of events with each multiplicity m as an independent 
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Poisson process. This yields a likelihood factor of the form 



•^neutron (^rad J '^muon) 

Pois(iVXk.' 



rad 



^muoii ) (1^) 



for each multiplicity m, where Pois (x | y) denotes the Pois- 
son probability to observe x events when y are expected, 
and N"^^ is the number of observed events with multi- 
phcity m. In our case, we have N^^^ — 2, N^^^^ — 1, and 
= for all other m > 1. 

Each number -/V^^ can be partly due to neutrons from 
a radioactive source and to muon-induced neutrons. We 
emphasize that, since the multiplicity spectra for the two 
types of neutron mechanisms are significantly different, 
a fit can distinguish the two contributions based on the 
observed multiplicities. It will ultimately choose the dis- 
tribution that fits best to the data. 

For neutrons from radioactive sources, the measure- 
ment of the coincidence rate has sufficient statistics 
(Fig.[7|) and we obtain the ratios r™^^ directly by dividing 
the contents of the respective bins. On the other hand, for 
muon-induced neutrons the statistics of coincident events 
is low and statistical uncertainties of the histogram in 
Fig. [T] need to be taken into account. We accomplish 
this by fitting the observed multiplicity spectrum with a 
(purely heuristic) exponentially decaying function. This 
function is then evaluated to obtain the ratios r^^^^- We 
directly add the likelihood of this fit to our total likelihood 
function and determine the parameters of the exponential 
simultaneously with all the other estimates. This way, 
their statistical uncertainty is automatically considered. 

While an exponential fit to the multiplicity spectrum 
is clearly not the most general form, it is evident that 
the abundance of coincidences decreases towards higher 
multiplicities. Our fit models this behavior and provides 
a reasonable approximation to the observed spectrum, in 
particular at the low multiplicities which are most relevant 
here. 



5.5 Lead Recoil Background 



quantities are characteristic of the implantation profile of 
a-emitters in the clamps and can thus be assumed to be 
universal if the underlying implantation mechanism is the 
same for all clamps. 

The unknown parameters of this background model are 
constrained by including the reference region in the lead 
band as introduced in Section WM The technical realiza- 
tion is identical to the one discussed for the a-background. 
Since the other known sources of events play only a neg- 
ligible role in this region, this is purely a reference mea- 
surement for the parameters of the lead recoil background. 
Moreover, the inclusion of the reference region again al- 
lows to treat the quenching factor of lead as a free fit 
parameter, which is automatically constrained by the ob- 
served reference events. The corresponding uncertainty is 
thus directly built into the likelihood function. 



5.6 WIMP Signal 

The density PxiE, y \ m^, cwn) of a possible signal due to 
coherent WIMP-nucleon scatterings is the sum of three 
components, one for each possible recoiling nucleus in 
CaW04. We calculate the expected recoil energy spec- 
trum for each component as a function of the WIMP 
parameters, using the usual standard assumption of an 
isothermal WIMP halo of density 0.3GeV/cm^, with a 
galactic escape velocity of Ucsc — 544km/s and an asymp- 
totic velocity of Voo — 220km/s. Effects of the nuclear 
substructure are taken into account by the Helm form 
factor as given in [11 , and the resulting energy spectra 
are finally convolved with the resolution of the phonon 
detectors to obtain the ultimately measured spectrum. 

In the expected light yield distribution, we take into 
account the uncertainties of the quenching factors for all 
three nuclei, approximating the errors given in Section [2. 3| 
by Gaussians. We note that the inclusion of these uncer- 
tainties has only a minor influence onto our results and 
even doubling the errors has no relevant effect. In the same 
way, varying the galactic escape velocity within its uncer- 
tainties is found to have negligibly small effects. 



In analogy to the a-background, we study the lead recoil 
background in a region where the Pb-band is free of other 
sources of events. As discussed above, this indicates that 
the energy spectrum of this background has a decreasing 
tail towards lower energies which we model by an expo- 
nential starting at 90 keV (the upper energy bound of the 
lead reference region), on top of a constant contribution: 



diV, 



Pb 



dE 



(E) 



Ad 

^Pb 



Pb 



exp 



E -90 keV 

decay ^ 



(13) 



In contrast to the simple fit of this function discussed 
in Section 14.41 we take into account the differences be- 



tween the detector modules here: The free parameter ^p^^ 
may be module-dependent and describes the different re- 
coil background rates in the individual detector modules, 
while we use the same spectral decay length -E^ocay ^^^^ 
constant background term Cpb for all modules. The latter 



6 Results and Discussion 

In this section, we summarize the results obtained from 
the maximum likelihood fit as introduced above. 



6.1 Resulting Fit Parameters 

We find that the total likelihood function has two max- 
ima in the parameter space, which we denote Ml and M2, 
respectively. Ml is the global maximum, but M2 is only 
slightly disfavored with respect to Ml. We will hence dis- 
cuss both solutions in the following. 

Table |4] shows the expected contributions of the back- 
grounds and of a possible WIMP signal in the two like- 
lihood maxima. The background contributions are very 
similar for Ml and M2: The expected e/7-background is 
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Ml 


M2 


6 /^-events 


8.00 ± 0.05 


8.00 ± 0.05 


Q-events 


11.5 tii 




neutron events 


7.5 It.? 


9.7 It} 


Pb recoils 


15.0tt.? 


18.7 tt? 


signal events 


29.4 i?;? 


24.2 1*;^ 


[GeV] 


25.3 


11.6 


awN [pb] 


1.6- 10"^ 


3.7- lO"'^ 



Table 4. Results of the ma:ximum likelihood fit. Shown are 
the expected total contributions from the backgrounds consid- 
ered as well as from a possible WIMP signal, for the parameter 
values of the two likelihood maxima. The small statistical er- 
ror given for the e/7-background reflects the large number of 
observed events in the e/7-band. The other errors correspond 
to a la confidence interval as determined by MINOS (see Sec- 
tion 5.1 1. The corresponding WIMP mass and interaction cross 



section are listed for each of the two likelihood maxima. 



one event per module according to the choice of the ac- 
ceptance region, vifith a negligible statistical uncertainty 
due to the large number of events in the e/7-band. The 
lead recoil and the a-background are similar to our simple 
estimates given in Section |4] Both these backgrounds are 
slightly larger than the contribution from neutron scatter- 
ings. In the context of the latter, the fit assigns roughly 
half of the coincident events to neutrons from a radioac- 
tive source and to muon-induced neutrons, respectively. 
This translates into about 10 % of the single neutron back- 
ground being muon-induced. 

In both likelihood maxima the largest contribution is 
assigned to a possible WIMP signal. The main difference 
between the two likelihood maxima concerns the best-fit 
WIMP mass and the corresponding cross section, with 

= 25.3 GeV in case of Ml and = 11.6 GeV for the 
case M2. The possibility of two different solutions for the 
WIMP mass can be understood as a consequence of the 
different nuclei present in our target material. The given 
shape of the observed energy spectrum can be explained 
by two sets of WIMP parameters: in the case of Ml, the 
WIMPs are heavy enough to detectably scatter off tung- 
sten nuclei (cp. Fig. [T]), about 69 % of the recoils are on 
tungsten, ~ 25 % on calcium and ^ 7 % on oxygen, while 
in M2, oxygen (52 %) and calcium recoils (48 %) constitute 
the observed signal and lead to a similar spectral distri- 
bution in terms of the recoil energy. The two possibilities 
can, in principle, be discriminated by the light yield dis- 
tribution of the signal events. However, at the low recoil 
energies in question, there is considerable overlap between 
the oxygen, calcium, and tungsten bands, so that we can 
currently not completely resolve the ambiguity. This may, 
however, change in a future run of the experiment. 

Fig. [11] illustrates the fit result, showing an energy 
spectrum of all accepted events together with the expected 
contributions of backgrounds and WIMP signal. The solid 
lines correspond to the likelihood maximum Ml, while 
the dashed lines belong to M2. The complicated shape 



total 

WIMP signal 
ybck 

Pb recoil bck 
a bck 

neutron bck 




25 

Energy [keV] 

Fig. 11. (Color online) Energy spectrum of the accepted 
events from all detector modules, together with the expected 
contributions from the considered backgrounds and a WIMP 
signal, as inferred from the likelihood fit. The solid and dashed 
lines correspond to the fit results Ml and M2, respectively. 



of the expectations is the result of taking into account 
the energy-dependent detector acceptances. In particular, 
the different energy thresholds of the individual detector 
modules lead to a steep increase of the expectations when 
an additional module sets in. 

We note that neither the expected a- or lead recoil 
backgrounds nor a possible neutron background resemble 
a WIMP signal in terms of the shape of their energy spec- 
trum. Even if our analysis severely underestimated one 
of these backgrounds, this could therefore hardly be the 
explanation of the observed event excess. 

On the other hand, the leakage of e/7-events rises 
steeply towards low energies and one may be tempted to 
consider a strongly underestimated e/7-background as the 
source of the observation. However, in addition to the en- 
ergy spectrum, also the distribution in the light yield pa- 
rameter needs to be taken into account. Fig. [12] shows the 
corresponding light yield spectrum of the accepted events, 
together with the expectations from all considered sources. 
Again, the shape of the expectations is the result of the 
individual detector acceptances being considered. As ex- 
pected, the contributions from the e/7- and also from the 
a-background quickly decrease towards lower light yields 
and thus differ significantly from the expected distribution 
of a WIMP signal. 

In order to check the quality of the likelihood fit, we 
calculate a p- value according to the procedure summarized 
in Section [O] We divide the energy-light yield plane into 
bins of 1 keV and 0.02, respectively, and include the accep- 
tance region of each module as well as the alpha- and Pb 
recoil reference regions in the calculation. The two likeli- 
hood maxima are found to give very similar results, with 
p- values of about 0.36 and 0.35, respectively. This not very 
small value for p indicates an acceptable description by our 
background-and-signal model. 
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light yield 

Fig. 12. (Color online) Light yield distribution of the accepted 
events, together with the expected contributions of the back- 
grounds and the possible signal. The solid and dashed lines 
correspond to the parameter values in Ml and M2, respec- 
tively. 

6.2 Significance of a Signal 




WIMP mass [GeV] 



Fig. 13. The WIMP parameter space compatible with the 
CRESST results discussed here, using the background model 
described in the text, together with the exclusion limits from 
CDMS-II [H], XENONIOO [l3j, and EDELWEISS-II [14], as 
well as the CRESST limit obtained in an earlier run [1]. Ad- 
ditionally, we show the 90% confidence regions favored by Co- 
GeNT [ig and DAM A/LIBRA tl6j (without and with ion 
channeling). The CRESST contours have been calculated with 
respect to the global likelihood maximum Ml. 



As described in Section [5T| the likelihood function can be 
used to infer whether our observation can be statistically 
explained by the assumed backgrounds alone. To this end, 
we employ the likelihood ratio test. The result of this test 
naturally depends on the best fit point in parameter space, 
and we thus perform the test for both likelihood maxima 
discussed above. The resulting statistical significances, at 
which we can reject the background-only hypothesis, are 



for Ml: 
for M2: 



4.7cr 
4.2 cr. 



In the light of this result it seems unlikely that the 
backgrounds which have been considered can explain the 
data, and an additional source of events is indicated. 
Dark Matter particles, in the form of coherently scatter- 
ing WIMPs, would be a source with suitable properties. 
We note, however, that the background contributions are 
still relatively large. A reduction of the overall background 
level will reduce remaining uncertainties in modeling these 
backgrounds and is planned for the next run of CRESST 
(see Section [?]). 



6.3 WIMP Parameter Space 

In spite of this uncertainty, it is interesting to study the 
WIMP parameter space which would be compatible with 
our observations. Fig. 13 shows the location of the two 



likelihood maxima in the (m^, CTwN)-plaiie, together with 
the 1(7 and 2a confidence regions derived as described in 
Section 15.11 The contours have been calculated with re- 
spect to the global likelihood maximum Ml. We note that 
the parameters compatible with our observation are con- 
sistent with the CRESST exclusion limit obtained in an 



earlier run pj, but in considerable tension with the limits 
pubhshed by the CDMS-II |12| and XENONIOO [I3] ex- 
periments. The parameter regions compatible with the ob- 
servation of DAMA/LIBRA (regions taken from Unj) and 
CoGeNT |15| are located somewhat outside the CRESST 
region. 



7 Future Developments 

Several detector improvements aimed at a reduction of the 
overall background level are currently being implemented. 
The most important one addresses the reduction of the al- 
pha and lead recoil backgrounds. The bronze clamps hold- 
ing the target crystal were identified as the source of these 
two types of backgrounds. They will be replaced by clamps 
with a substantially lower level of contamination. A sig- 
nificant reduction of this background would evidently re- 
duce the overall uncertainties of our background models 
and allow for a much more reliable identification of the 
properties of a possible signal. 

Another modification addresses the neutron back- 
ground. An additional layer of polyethylene shielding 
(PE), installed inside the vacuum can of the cryostat, will 
complement the present neutron PE shielding which is 
located outside the lead and copper shieldings. 

The last background discussed in this work is the leak- 
age from the e/7-band. Most of these background events 
are due to internal contaminations of the target crystals 
so that the search for alternative, cleaner materials and/or 
production procedures is of high importance. The mate- 
rial ZnW04, already tested in this run, is a promising 
candidate in this respect. 
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8 Summary 

With an exposure of 730 kg d, the CRESST Dark Mat- 
ter search has observed in the latest run a total of 67 
events in an acceptance region of low energy nuclear re- 
coils. Possible background contributions to this number 
include leakage of e/7-events, events from neutrons, from 
alpha-particles, and from recoiling nuclei in a-decays. We 
have estimated these four backgrounds and have found 
using a likelihood ratio test that, at a significance larger 
than 4cr, these backgrounds are not sufficient to explain 
all the observed events. Scatterings of WIMPs may be the 
origin of this effect and, under this assumption, we have 
derived the corresponding WIMP parameters. Finally, we 
have presented the plans for the next run of the experi- 
ment, in which we aim for a further clarification of this 
hypothesis. 
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